Trapping and Wiggling: Elastohydrodynamics of Driven Microfilaments 
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We present a general theoretical analysis of semiflexible filaments subject to viscous drag or point 
forcing. These are the relevant forces in dynamic experiments designed to measure biopolymer 
bending moduli. By analogy with the "Stokes problems" in hydrodynamics (fluid motion induced 
by that of a wall bounding a viscous fluid), we consider the motion of a polymer one end of which 
is moved in an impulsive or oscillatory way. Analytical solutions for the time-dependent shapes of 
such moving polymers are obtained within an analysis applicable to small-amplitude deformations. 
In the case of oscillatory driving, particular attention is paid to a characteristic length determined 
by the frequency of oscillation, the polymer persistence length, and the viscous drag coefficient. 
Experiments on actin filaments manipulated with optical traps confirm the scaling law predicted by 
the analysis and provide a new technique for measuring the elastic bending modulus. A re-analysis 
of several published experiments on microtubules is also presented. 



I. INTRODUCTION 

The attempts of theoretical physicists to contribute in 
some useful way to the study of biology has, so far, been 
most successful in systems in which all forces and motion 
can be modeled and mathematized explicitly, or those 
governed by equilibrium statistical mechanics, for which 
equipartition can be invoked. One specific example of 
this success is the analysis of structural microfilaments, 
essentially one-dimensional mechanical objects with no 
moving parts. Despite this unassuming mechanical de- 
scription, these semiflexible biopolymers are responsible 
for innumerable essential functions at the molecular and 
cellular level. 

Depending on the bending modulus of the filament 
in question, experiments on semiflexible biopolymers 
largely rely on either mechanical or statistical techniques. 
Microtubules, with a persistence length ~ 5 mm, are 
quite amenable to micromanipulation or forcing via hy- 
drostatic drag. Actin and nucleic acids, with persistence 
lengths near and 50 nm, respectively, fall in the realm 
of statistical mechanics (note that we are here address- 
ing the bending elasticity, not the stretching elasticity, 
another area of great excitement and successes [5, 30, 
37]). 

However, much of the analysis of elasticity in this con- 
text seems to ignore a fact known to the ancients who 
developed the field: boundaries matter for finite objects. 
While often in physics boundary effects can be ignored, 
the fact that biopolymers are finite has important and 
unexpected consequences and must be included carefully 
the analysis. These effects appear because dynamics of 



the elastica* are governed by fourth-order spatial deriva- 
tives. This fact forces us to rethink what are the most 
appropriate basis functions in which to model dynam- 
ics. While this question may seem purely academic, we 
will see that it is quite crucial in attempts to develop 
accurate dynamic methods for the quantitative study of 
semiflexible biopolymers 

Speciflcally, we couple elasticity theory and over- 
damped viscous hydrodynamics (as is appropriate in 
the biological context) to explore elastohydrodynamics. 
While experiments coupling elasticity and drag have been 
presented before [10, 33], treatments of dynamics, rather 
than mechanical equilibrium, have been less than com- 
plete. Although equations with the appropriate units will 
be sufficient to determine the scales of forces and veloci- 
ties, if we wish to extract numbers from the experiments 
it is necessary to have a thorough analysis. The slen- 
derness of the filaments allows us greatly to simplify the 
hydrodynamics and arrive at a local partial differential 
equation of motion. Moreover, we will find that coupling 
to hydrodynamics allows us to extend the range of me- 
chanical experiments to much smaller bending moduli. 
For example, whereas measurements of actin's rigidity so 
far have been via fluctuation analyses invoking equipar- 
tition and thus statistical mechanics, we present here an 
experimental method which does not rely on nonzero 
temperature. Further, the method allows investigation 
into questions which have been raised as to whether or 
not actin can even be treated as a semiflexible polymer, 
or is in fact scale-sensitive [16] or dynamic in its elastic- 



*a line of particles for which the resistance to bending is a 
couple proportional to the curvature [21] 
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ity. Such a purely mechanical treatment obviates pos- 
sible complications to statistical treatments like dimen- 
sionality [26] , correlations among sampled images, or self- 
avoidance. 

It is our hope that this new experimental method, as 
well as the general analytic techniques here outlined, will 
contribute to the current exciting and active dialogue 
among physicists and biologists regarding the nature and 
numbers behind biopolymcr rigidity as well as the effects 
of varying proteins and environments on elastic moduli. 
Further, these new methods and analyses should prove 
useful in the study of other examples of dynamic elastic 
filaments, such as supercoiled fibers of B. subtilis [24]. 
We intend this investigation to be the necessary starting 
point for such exciting extensions. 
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velocity at any point x and time t depends only on the 
ratio x/{utY/'^; likewise, in the oscillatory case, oscilla- 
tions decay with a characteristic length that scales as 

We introduce here the analogous two problems in elas- 
tohydrodynamics, illustrated in Fig. 2. They involve (I) 
the deflection of a polymer anchored at one end following 
the instantaneous introduction of a uniform fluid veloc- 
ity U , and (II) the steady undulations of a polymer one 
end of which is oscillated. Rather than a diffusion equa- 
tion as in the Stokes problems, the dynamics of small 
deformations y{x,t) of the filament are governed by a 
fourth-order PDE of the form 

Vt = -i>yxxxx , (1-2) 

where i) = A/( plays the role of a "hyperdiffusion" coeffi- 
cient, with A the bending modulus and ( the drag coeffi- 
cient. This equation has appeared before in the literature 
on semiflexible biopolymers [1, 12], primarily in the con- 
text of scaling arguments for relaxation times; our goal 
here is to provide a complete solution given arbitrary ini- 
tial and boundary conditions as dictated by experiment.^ 
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FIG. 1. Geometry of Stokes problems I and II. 

A useful starting point for developing the general 
mathematical treatment of the dynamics of an elastic 
filament in a viscous medium will be to consider the sim- 
plest time dependencies possible. To that end, recall the 
classic problems introduced by G.G. Stokes in 1851 [31], 
illustrated in Fig. 1, involving the motion of a viscous 
fluid bounded by a wall that is either (I) moved impul- 
sively or (II) oscillated. These easily-solvable problems 
capture the essential ideas of viscous diffusion of veloc- 
ity. The experimental geometry is such that the Navier- 
Stokes equation for the velocity field u[x, t) is simply the 
diffusion equation 

Ut = VUxx , (1-1) 

where v — ^1 pis the kinematic viscosity, and /i and p are 
the fiuid viscosity and density. Subscripts on functions 
will indicate differentiation throughout unless otherwise 
indicated. The salient features of the solutions are the 
relationships between length scales, time scales, and ma- 
terial parameters. Specifically, in the impulsive case, the 
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FIG. 2. Geometry of elastohydrodynamic problems I and II. 

An analysis similar to that presented below of the oscil- 
latory passive elastica was carried out a number of years 
ago by K. E. Machin [22, 23], who considered the motion 



' Nota Bene: in [1], Eq. 2 should have a minus sign before 
the bending constant; as written, the equation is ill-posed. 
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of a driven flagcllum. Machin was interested specifically 
in a semiinfinite active flagellum which was bent with a 
set of boundary conditions amenable to analysis. Ours 
will be more malicious, but not subtle. 

We will first derive the general equations of motion for 
elastica embedded in viscous flow. We then apply this 
dynamic to a number of scenarios. Inspired by "Stokes 
problems I and 11" in fluid dynamics (SI and SII), we first 
solve "problems I and 11" of elastohydrodynamics (EHDI 
and EHDII), each of whose dynamic mimics its hydrody- 
namic analogue. Problem I will require a bit of math- 
ematical construction to assist our physical intuitions. 
Specifically, we must construct a set of basis functions ap- 
propriate to the equation of motion and specified bound- 
ary conditions. All the pleasant features found when ap- 
plying Fourier space to unbounded or periodic systems 
are found here as well, in what we term W-space. Unlike 
Fourier space, W-space respects both the compact sup- 
port and boundary conditions of the elastica and thus 
diagonalizes the equation of motion. We then discuss 
an experimental realization of problem II and its analy- 
sis, which provides a new technique for the measurement 
of a polymer's bending modulus. Finally, we comment 
on some recent related experimental work by a separate 
group, brought to our attention as the original version of 
this paper was being completed. 



II. ELASTIC FORCES 

A bent elastic polymer exerts a restorative force given 
by the functional derivative of a bending energy 



which for an elastic curve may be expressed as' 
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(2.1) 



(2.2) 



Here, A is the bending stiffness constant, with units of 
energy x length. It may also be expressed as the prod- 
uct EI of the Young's modulus E and the moment of 
inertia / [11]. For a polymer of persistence length Lp at 
absolute temperature T, exploring all configurations in 
D dimensions, we may also derive by equipartition the 
equivalence A = {D — l)kBTLp/2. 

Henceforth we will consider elastic filaments that lie in 
a plane. This is the geometry most well-suited to data- 
taking via microscopy. The curvature n may then be 



expressed exactly as dd/ds, where 9 is the angle between 
the tangent to the curve and some fixed axis (see Fig. 
3); the .x-axis is a convenient choice, for which tan^ = 
{dy / ds) / {dx / ds) = dy/dx, which implies 



K = 



3/2 



(2.3) 



Taking the functional derivative of the energy (2.2), we 
derive the force and boundary conditions 



fs=A ^Kss + ^K^j n , 
K = Ks = at free ends. 



(2.4) 
(2.5) 



These boundary conditions indicate the torquclcssness 
and forcelessness, respectively, at free ends of an elas- 
tica [34]. At hinged or clamped ends different boundary 
conditions hold, as will be discussed below. 
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FIG. 3. Geometry of an elastic filament: il=local radius of 
curvature; t, n=unit tangent and normal, 9 = cos~^(t • ex); 

d=diametcr of the filament; arclength s varies from to L, the 
total arclength. Within the approximations of slender-body 
hydrodynamics, a local anisotropic proportionality is satisfied 
between an external force per unit length f and the velocity 

V. 

We may linearize these results for small deviations 
from a horizontal line, {yx <C 1) or equivalently linearize 
Eq. 2.2 to find k ~ y^x + O (yl) and 



fe -Ayxxxxey + 0{yl) , 



(2.6) 



(where ey is the unit vector in the y direction) with 
boundary conditions 



yxx = and y^ 



at free ends. 



(2.7) 



The specification of the filament dynamics is complete 
upon definition of the hydrodynamic drag which balances 
ff. We now turn to this problem. 



■'■Note that we may also include any forces of constraint, 
such as a Lagrangian tension to enforce inextensibility [13], 
but such terms will be of higher order in the curvature than 
we will consider in this investigation. 



III. SLENDER-BODY HYDRODYNAMICS 

We wish to consider experiments taking place on cel- 
lular biological scales, with typical lengths L in microns. 
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fd = C [nfi + /3tt] • (r^ - u) . 



times t in seconds, and dynamic viscosity /i that of water, 
in centipoise. Tlie nondimcnsional parameter character- 
izing the hydrodynamic behavior, the Reynolds number, 
is [/L/i/ ~ L'^/tf ~ 10-710-2 ~ 10-6. We are thus 
safely in the low-Reynolds number or 'Stokesian' regime. 
In this Aristotelian overdamped limit, forces balance ve- 
locities rather than accelerations. For a body whose 
length is much greater than its width, the well-developed 
set of calculations known as 'slender-body hydrodynam- 
ics' applies [6, 7, 17]. If this filamentous polymer has 
diameter d, length L, and an aspect ratio rf/L <C 1, we 
have the simplified, local, anisotropic proportionality be- 
tween the drag force and the velocity r^, 

(3.1) 

Here, (s) is a force per unit length exerted on the fila- 
ment, and n (s) and t (s) are unit vectors in the normal 
and tangential directions at arclength s along the poly- 
mer. The velocity of the polymer is denoted rj (s), and u 
is any background velocity which may be present in the 
problem; the drag should be a function of the former rela- 
tive to the latter. The anisotropy, evident when dragging 
a pencil through molasses, between motions parallel and 
perpendicular to a slender object's long axis is embodied 
by the parameter /3, which depends logarithmically on 
the aspect ratio, with asymptotic behavior 

lim /3 = - . (3.2) 

An active field of applied mathematics in the last twenty 

years has been calciilations of the viscous drag coefficient 
( for various slender shapes [4, 6, 7, 17, 20, 29]. The 
limiting behavior for small d/L has the form 



c = 



In iL/d) + c 



(3.3) 



where c is a constant of order unity which depends on 
the shape of the body. 

Now that we have defined the hydrodynamic drag, we 
demand that it balances elastic forcing (f^ = f^) to derive 
the equation of motion: 



C [nn + pit] • (rt - u) = A 



1 . 



(3.4) 



Linearizing the expression for the drag (3.1) for nearly- 
straight polymers {yx «C 1), and noting that its tangen- 
tial components are of order y^, the dynamic reduces to 



C {yt -u) = -Ayx 



(3.5) 



Here, u = u-ey. In the absence of any background flow, 
Vt = -^Vxxxx, (3.6) 

where D = A/C, as \n Eq. 1.2. This is the simplest lin- 
earized expression of elastohydrodynamics: elastic forces, 
characterized by a fourth spatial derivative, balance vis- 
cous drag. It shares many similarities with the diffusion 
equation (1.1), and may be thought of as "hyperdiffu- 
sion" of displacement in analogy with hyperviscosity. 



IV. ELASTOHYDRODYNAMIC PROBLEM I 

Now that we have established the equation of motion 
appropriate to these elastohydrodynamic analogs, we re- 
call the solutions to the fluid dynamics problems SI and 
SII in hopes of exploiting the analogy as much as pos- 
sible. In Stokes I (SI), a semiinfinite plane of fiuid is 
forced by a wall which is motionless for time t < and 
has velocity Uey for t > 0. In Stokes II (SII), the wall 
oscillates as U cos {ujt) ey, and we solve for the behavior 
after transients have died away. 

As illustrated in Fig. 1, velocity gradients are in the 
a;-direction in both cases, and hence perpendicular to the 
direction of flow (along the y-axis). In the absence of an 
imposed pressure gradient, the Navier-Stokes criuation 
for the fluid velocity u (x, t) parallel to the wall is simply 
the linear diffusion equation (1.1), Ut = yuxx- 

A convenient method of solving SI with the associated 
boundary condition is to postulate a scaling solution in- 
spired by dimensional analysis: u{x,t) = UF{^), with 
1/2 

^ = x/ (vt) ' . The scaling function F then obeys a 
nonautonomous ordinary differential equation 



(4.1) 



The solution to this \s F = erfc (^ / \/2) , where crfc is 
the complementary error function. The dedicated may 
arrive at the same result more laboriously via Laplace 
transforming Eq. 1.1. 

Armed with some understanding of SI, we now turn to 
problem I of elastohydrodynamics (EHDI). In problem 
I, we consider an elastic filamentous polymer which is 
anchored at the origin. For i < it lies along the line 
segment {y = 0; < a; < L}. We then may consider 
forcing the filament by moving one end relative to the 
fiuid (moving the anchor) or moving the fiuid relative to 
the polymer (moving, for example, the cover slip). We 
will first attempt to do this in a way as analogous to SI 
as possible. 



A. Semiinfinite polymer 

Consider first moving the anchor relative to the fluid 
with velocity u = uey. This may be thought of as mov- 
ing a bead optically trapped to one end, or some other 
micromanipulation. We now observe that the equation 
of motion, yt = —vyxxxx, inspires a scaling solution of 
the form y = utF{^), where ^ = x/ [vt)^^^. Inserting 
this special form of the solution, the equation becomes 



(4.2) 



an equation not unlike Eq. 4.1. 

The solution F must admit the conditions that it and 
its derivatives approach for a fixed time as a; — > 00, so 



4 



lini{F(0,Fj(0,%(0,---}=0. (4.3) 

Moreover, specifying the anchor position impHes 

ut = y (0, t) = utF (0) F{0) = 1 . (4.4) 

Depending on whether we wish to describe a filament 
anchored by, for example, an axoneme or a bead, we 
demand (0) = or y^x (0) = 0, implying Fj(0) = or 
F^^ (0) = 0, respectively. 
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FIG. 4. Scaling solution to (1.1) for a semiinfinite clamped 
polymer. The meaning of "scaling solution" is that y, x, and 
t, are suitably rescaled such that, if a friend were to imprint 
this image on silly putty, and run towards you with a velocity 
~ t"®/** while stretching the putty in the y direction a length 
you would see the predicted shape of the semiinfinite 
elastica. For the most accurate results, have your friend start 
infinitely far away. This complication is obviated by consid- 
ering the finite elastica. 

Note that, had we instead chosen to consider a fixed 
anchor and a moving fluid, wc would merely implement 
two changes. First, the equation of motion becomes 
{yt — u) = —yyxxxx- Second, the boundary condition be- 
comes y (0) = 0. However, introducing a Galilean change 
of frame y' {x,t) = y — ut, we see that wc recover the 
equivalent equation of motion, boundary conditions, and 
thus problem. This relationship between boundary con- 
ditions and equations of motion will allow us to apply 
one analysis to a number of experiments. We will use the 
terms "homogeneous" and "inhomogeneous" referring to 
both boundary conditions and equations of motions since 
we may transform from one to the other so simply. 

Though nonautonomous, Eq. (4.2) is linear in F and 
quite amenable to numerical solution. A plot of the so- 
lution for the clamped (?/x(0) = 0) polymer is shown in 
Fig. 4. Clearly, the function satisfies our qualitative intu- 
itions about the shape: an initially straight curve slopes 



towards as the argument grows, where all derivatives 
vanish. The appearance of the wiggle leads us more care- 
fully to consider the realm of validity of our solution. 
Given our scaling ansatz y = utF{^), we see that the 
slope yx = utF^/{v'ty^^ grows in time without bound, 
thus failing to meet the criterion on which the lineariza- 
tion of Eq. 3.4 was predicated: yx 1. Intuitively, we 
suspect that the neglected nonlinear terms would sup- 
press these undulations. 

Noting that \F^\ attains its largest value of ~ 0.51 for 
^* ~ 1.1, we see that we can only trust the solution to 
initial times such that 



,.1/3 



iuF^C)) 



4/3 



2.5 



1/3 



(4.5) 



For an experimentally convenient velocity of ~ 10/x/s, 
this time is order a fraction of a second for actin, and 
a second for microtubules. Wc have thus seen that the 
obvious extension of the hydrodynamic analogue is not 
necessarily the correct approach - a conclusion which will 
be reached repeatedly in this investigation. 

Enjoyable though it is to contemplate the long-time 
behavior of the semiinfinite dragged elastica, a more prof- 
itable inquiry would be the experimentally realizable case 
of a finite polymer. We therefore turn now to the case 
of the finite elastica, clamped at one end and free at the 
other, and subject to impulsive hydrodynamic drag. Our 
analysis is applicable to experiments in which either the 
cover slip is moved or, as a special case, in which the 
elastica is allowed to relax from some initial conditions 
in the absence of flow. 



B. Finite polymer 

In order to make the mathematics as transparent as 
possible, wc first nondimcnsionalizc the equation of mo- 
tion (3.5). Distances in x are rescaled by the total length 
L, time by the elastohydrodynamic time scale (L'^/A, 
and distances in y by the velocity times this time scale: 



X = aL , t = T 



y{x,t) 



u — —h (a, r) 



(4.6) 



The governing equation, [y^ — 
comes 

(hr - 1) = -ho 



u) 



-^Vx 



then be- 



(4.7) 



Like any other inhomogeneous differential equation, the 
general solution of (4.7) must be the homogeneous solu- 
tion {hr = —haaaa) plus thc particular solution (/i-r — 1 = 
—haaaa)- We indicate this by defining 



where g and / satisfy 



9r 



-9a 



f + g 



It 1 fa.c 



(4., 



(4.9) 



We will find that these two functions have markedly dif- 
ferent qualitative behavior. 
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1. Homogeneous solution 

The homogeneous equation is Qt = —Qaaaa- Well- 
vcrscd in the litany of Fourier transforms, wc first left- 
multiply by an as-yet arbitrary function Wu (a) (where 
k indicates a parameter rather than a derivative) and 
integrate over the domain of a, 



daWkdrg = - / daWkd^g . 
'0 Jo 

Integrating by parts we obtain 



(4.10) 



[ daWk9 = - I dadt{Wk)g 
Jo Jo 

- mdlg\l+diWkdlg\l 

- dlWkd'^g\l+d'^Wkg\l. 



(4.11) 



Observe that the integration by parts of the fourth- 
order derivative introduces 8 separate surface terms. The 
boundary conditions implied by the functional derivative 
(2.5) dictate the vanishing of the second and third deriva- 
tives at the free end {x = L). Requiring g to satisfy this 
condition eliminates two of the 8 terms. 

The left end of the polymer is clamped at the origin, 
so y{x = 0) = yx {x = 0) = 0. Demanding this behavior 
of g eliminates two additional surface terms. We now 
choose yVk to satisfy the same boundary conditions as y 
and g, and therefore h: 



Wfe(0) = 0, 9„Wfe(0)=0, 

dlWk (1) = , a^Wfc (1) = 0. 



(4.12) 



This annihilates the remaining 4 surface terms in (4.11). 
Finally, we choose Wfe to obey 



(4.13) 



Defining gk = J^daWkg , the equation of motion be- 
comes d^gk = —k^gk, the solution to which is 



gk (t) = gk (0) e" 



-fc*r 



(4.14) 



If we wish to describe the dynamics in such terms, we 
must construct the Wk, which necessitates that we iden- 
tify the allowed values of k. We turn now to this problem. 



as well in the basis of the function space. The general 
solution of Eq. 4.13 is 



Wfe = ai sin [ka) +a2 cos (fca) 
+03 sinh {ka) +04 cosh {ka) . 



(4.15) 



The expression has four unknowns, as a solution to a 
fourth-order problem must. The a = (left end) bound- 
ary conditions yield 04 = —a^ and 03 = — ai. The last 
two conditions, d^Wk (1) = and d%Wk (1) = 0, can be 
expressed as a matrix equation 



/ sinh A; -h sin fc cosh k + cos k 
\ cosh k + cos k sinh fc — sin fc 



0. (4.16) 



Since the matrix has a zero eigenvalue, it must be singu- 
lar. This solvability condition is 

1 



cosfc 



cosh A; 



(4.17) 
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FIG. 5. Solution of the transcendental equation for the 
roots fcn- (a) Graphical construction of the solution, (b) The 
roots (circles) asymptote to fcn = tt (2n — 1) /2 (dashed line). 



2. Explicit construction of the Wk 

A moment's thought reveals that the Wk can not be 
constructed out of simply the familiar sin's and cos's 
of Fourier space, which are incompatible with boundary 
conditions in which successive derivatives vanish (c/. Eq. 
2.5). A countably infinite family of such Wk can, how- 
ever, be constructed by including hyperbolic functions 



As Fig. 5 illustrates, this transcendental equation has 
an infinite number of solutions. For large values of k, as 
1/ coshfc 0, the solutions approach the solutions of the 
Fourier-like solvability condition cos/c = 0, i.e., kn+i — » 
^ -I- irn. The first few solutions are 

37r 

A;i ~ f 4- .304 ~ 1.875 ,k2^ — ~ .018 ~ 4.694 , 

fcs ~ -I- .001 ~ 7.855, fc4 ^ Y - ^^-^^^ ' • • • (^-^^^ 
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We can insert either row of the LHS of (4.16) to solve 
for the relationship between 02 and ai, which we will 
express here as 



02 = - 



sin k + sinh k 
cos k + cosh k 



ai. 



(4.19) 



The constants themselves are functions of k. Note also 

that as A: ^ ±00, 02 — > ±1. We now choose ai such that 
Wfe is normalized, i.e., daW^ = 1. The first three 
eigenfunctions are shown in Fig. 6. 

Wc have now completely and explicitly constructed the 
family of functions Wk which are normalized eigenfunc- 
tions of with particular boundary conditions. Note 
that had we chosen some other set of boundary condi- 
tions, a different solvability condition and eigenfamily 
would have resulted (c/. Appendix A). For example, 
in the case of the elastica with free ends, we employ an 
expansion of y with the basis functions of Eq. (A6). 
Note that if we differentiate this expression and define 
C(n = kn/2, then we recover Eq. 28 of [12]. The expan- 
sion in terms of k has the advantage of a single solvabil- 
ity condition cos(A'„) cosh(fc„) = 1, rather than the two 
conditions tanh(Q;„) — (— 1)" tan(a„). The latter two 
conditions may be derived from the former via half-angle 
formulae. 

This operator is self-adjoint (c/. Appendix B) and thus 
the eigenfunctions constitute a complete basis in function 
space onto which we may project initial data and relate to 
later-time solutions via Eq. 4.14 in the standard Green's 
function way: 



g{a,T) =J2\^k)m\g{a,r)) 

k 

fe 

= / da'Q {a,a']T) g{a' ,Q) , 
Jo 

where the Green's function is 



(4.20) 



(4.21) 



This is the exact solution of the linearized homogeneous 
equation. Note that the compact support and the bound- 
ary conditions break translation invariance, reflected in 
the fact g can not be expressed as Q {a — a'^r). 
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FIG. 6. The first three eigenfunctions for BHD problem I. 
The dotted line indicates the normalized third-order polyno- 
mial describing an elastica bent by a point force at the right 
end. Note the surprising overlap with Wfci, as will be ex- 
ploited in section 4. 

We note from the exact solution (4.21) that each mode 
gk decays independently and exponentially with time. 
This is to be compared with diffusive problems, in which 
each mode decays exponentially in time, except for the 
zero (average) mode, which is constant. In this experi- 
ment, the boundary conditions arc incompatible with the 
existence of a zero mode. The system "hyperdiffuses" to 
homogeneity. 



3. Particular solution 



Since g{a, r) decays to zero as t c», we have 

lim h{a,T) =/(a,oo). (4.22) 

r— >oo 

Thus the particular solution, if it exists, must be the long- 
time solution; no matter what initial data, the solution 
to the inhomogeneous equation of motion (4.7) must be 
the global attractor. 

By the same argument, we must only construct a so- 
lution for it to be the correct one. We first look for a 
time-independent hmction, for which Eq. 4.9 becomes 
1 = faaaa- The Unique solution consistent with the 
boundary conditions (4.12) is 



(4.23) 
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4- Digression on boundary conditions 

It is useful to pause momentarily and reflect on the im- 
portance of boundary conditions and how they will enter 
into the analysis. To this end, we consider an analogous 
experiment in which the left side of the polymer is hinged 
rather than clamped. This would be realized by holding 
the polymer in an optical trap, for example, rather than 
some torquc-cxcrting anchor like an axoneme fixed to a 
cover slip. For this experimental setup we replace the 
boundary condition ha (0, r) = with haa (0, t) = 0. 

Intuitively, we expect the polymer to align itself with 
the flow as t ^ 00 in the absence of any torque at the 
left end. Mathematically, we may think of the change in 
boundary conditions via some curious pathologies. The 
flrst complication is that there does not exist a time inde- 
pendent fourth-order polynomial in a which is consistent 
with the boundary conditions. This prevents us from 
constructing a static attractor for the problem. How- 
ever, we note that the equation of motion is solved by a 
flfth-order, time-dependent polynomial in a: 
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(4.24) 



To derive this poynomial, we flrst express the solution 
in the form c„ (r) a". However, we are considering a 
Stokesian dynamic, in which the time dependence is first- 
order and the equation of motion is linear; the driving is a 
constant in time and thus time should enter only linearly 
into any steady state solution for the position (c/. [19]). 
Constraining dfcn = and respecting the relationship 
between 9^c„ and c„+4 dictated by Eq. 4.7, we uniquely 
specify the polynomial up to the constant C. 

The fact that our long-time solution has an arbitrary 
constant should lead us to rethink splitting the solution 
into a polynomial and a set of only-decaying modes. Re- 
turning to the set of eigenvalues of d^, we discover a sec- 
ond complication for this new boundary condition: there 
now exists a zero mode - a nontrivial solution, consistent 
with the boundary conditions, to the equation d^Wo = 0, 
i.e., the normalized polynomial ^/Sa. We then may 
choose C to eliminate the overlap of Wo with the initial 
data, i.e., = go = {g\Wo} = {h{a,0)\Wo} - {Ha)\Wo)- 

We now sec how the change in boundary conditions 
creates drastically different physical behavior. The long- 
time solution contains a term describing a straight line 
whose slope grows with velocity 3u/2L without bound. 

In order further to illustrate the relationship between 
a change in boundary conditions and the qualitative be- 
havior, note that the dynamics, even in the presence of an 
inhomogeneous equation of motion, can be cast in terms 
of the functional derivative of an energy: 



hr 



' 6h (a) ' 



[ da' 
Jo 



1 2 



(4.25) 



The first term represents the drag force acting in the 
positive y direction, while the second is simply the nondi- 

mcnsionalization of the elastic bending energy term from 
which we originally derived the equation of motion. We 
then find 



f dahrj^ =- [\a {Kf , (4.26) 
Jq oh (a) Jq 



indicating .F is a monotonically decreasing function in 
time [8, 13]. 

We must now only show that for the clamped (hinged) 
polymer this energy functional is (is not) bounded from 
below. Rewriting h = X]c„(r)a", we can evaluate the 
energy explicitly as 



71=0 



n + l 



+ n 



where 



n{n — 1) m (m — 1) 



m,n=2 



(4.27) 



(4.28) 



If h (0) = /i^(0) = 0, Co = ci = 0. Eq. 4.27 then becomes 



E 

n=2 



n+l 



(4.29) 



We see that we can not simply make the first term arbi- 
trarily negative by introducing a large and positive Cj for 
some j, since this term will appear quadratically (and al- 
ways positively) in the second term. However, replacing 
the condition ha (0) = with haa (0) = changes the 
condition ci = to C2 = 0, and the energy becomes 



n— J 



+ n . 



(4.30) 



Now the energy can become arbitrarily negative if Ci = 
h'{0) becomes arbitrarily positive. A divergent slope sim- 
ply means the curve points straight up, in accord with 
our intuition for the long time behavior of a polymer free 
to rotate in some background flow. Note such a long-time 
behavior means leaving the small-/ix limit for which the 
dynamic was originally derived. 



5. General solution 

Returning to the clamped polymer in the presence of 
some background flow, we project the definitional state- 
ment h {a, t) = g (a, t) +h (a) onto the Wfe (a): 



hk (t) = gk (t) -I- hk, 



(4.31) 
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which imphes the initial condition (0) = hk (0) — hk- 
Recalhng the simple time-dependence of the modes gk 
from Eq. (4.14), we see 



hk (r) = hk{l- e"'^'") + hk (0) e" 



(4.32) 



The dynamic thus mimics that of a capacitor, charging 
up with the final shape-state and draining of the initial 
shape-state, each mode governed independently with de- 
cay rate k~'^. In real space, 



h (a, r) =h{a)+Y^ Wk (a) e" 



fe*r 



k=ki 



X / da'Wk {a') [h {a', 0) - h {a')] . (4.33) 
Jo 

In the experiment considered, the initial condition is a 
flat polymer: h{a,T = 0) = 0. Since h is the solution to 
haaaa = 1, with boundary conditions h (0) = ha (0) = 
haa (1) = haaa (1) = 0, we find upon integrating by parts 
that 



^ da'Wk {a') h {a') = ^ 



(4.34) 



where Wk = Jq daWk, and thus Eq. 4.33 reduces to 



h (a, T) = h{a)-f^ Wk (a) ^e''^'^ 

k=ki 



(4.35) 



Evaluating the first few integrals, we find for hk = 
J daWkh = Wk/k\ 

hk,^6x 10-2 , /ife, ~ 9 X 10-^ , hk,^7x 10-^ , 
hk,^lx 10"^ . . . (4.36) 

Each mode with k > ki decays exponentially faster in 
time than the lowest mode, which thus dominates as r ^ 
oo, so 



h h{a) - hk.e-"^^ ~ /i - QMWk.e''^-''''^ ■ (4.37) 

Our picture of the impulsive dynamic of elastica in 
viscous flow is thus as follows: wc project onto a special 
function space in which the long-time solution and the 
difference between initial data and the long-time solution 
exponentially charge and decay, respectively, each mode 
behaving independently. We arc left with only the the 
long-time solution as the asymptotic limit r ^ oo. 



V. ELASTOHYDRODYNAMIC PROBLEM II 

In Stokes II, the driving force is exerted by a wall os- 
cillating with velocity u = C/ey cos (wt), or position y = 



Do cos (ojt) . To solve the steady-state limit of SII, we pos- 
tulate u{x,t) = C/K(e*'^*G (??)), where rj = x{u)lvf''^. 
Inserting into Eq. 1.1, we then see G satisfies 



(5.1) 



for which the solution vanishing as ?7 ^ oo is G = e 
We then find u {x,t) = Ue~'^^^ cos [ut — rj/V^), or, in 
a form useful for comparison to the elastohydrodynamic 
case, 



u {x, t) = Ue-^*'' cos ((747? - tot) 



(5.2) 



where C4 = cos(7r/4) and ^4 = sin(7r/4). This solution 
describes right-moving waves of velocity uii/C^, decaying 
as X ^ 00 with decay length 1/84. 

We now consider a polymer held by an optical trap 
which moves with position y [x = 0) = yo cos {ivt). Since 
we have shown in the previous section that all modes 
satisfying the homogeneous equation of motion with ho- 
mogeneous boundary conditions decay exponentially, we 
must only find a solution in the presence of inhomogene- 
ity (here, the driving) to find the long-time limit of the 
dynamic. 

In order to verify the validity of our analysis as well 
as the plausibility of EHDII as a method for measuring 
biopolymcr rigidity, we conducted the experiment and 
analyzed image data as described below. A scaling re- 
lation predicted by the analysis was confirmed, and a 
new method for measurement of the persistence length 
of actin was demonstrated. 

The experimental setup is shown in Fig. 7: F-actin 
is bound to a silica bead which is optically trapped. By 
oscillating the position of the bead sinusoidally in time, 
the filament wiggles back and forth, propagating waves of 
displacement down its length. The motion relative to the 
fluid is opposed by the fluid viscosity, and the "wiggles" 
are opposed by the elasticity of the polymer. 



slide 



bead 



coverslip 



20(J. 



oil 



• microscope 
/ objective 



Nd: 


YAG laser beam 




>^1.054|J, 



FIG. 7. Experimental setup. 



9 



A. The chciracteristic length 

The elastic constant A has units of energy x length, 

while the viscoiis force per unit length per unit velocity 
has the dimensions of a viscosity or action density /z: 



[C] = [m] 



mass 



length X time 



energy x time 
length^ 



(5.3) 



Thus, the natural length obtained from A, (, and the 
frequency of oscillation u is 



A 



1/4 



1/4 



(5.4) 



Nota bene that is not a mere rescaling of the persis- 
tence length. 

With a previously-published persistence length for 
actin of Lp ~ 15/Ltm [26], a viscosity fi = 0.01 cp, 



ksT ~ 4 X 10"" erg at T 
CO in units of sec~^, we obtain 



300 K, and measuring 



(5.5) 



Thus for frequencies on the order of 1 Hz, we obtain 
length scales of order microns, somewhat below the per- 
sistence length. This range of frequencies seems quite 
advantageous for experiment. 

This elastohydrodynamic length £ (lu) is precisely the 
length found upon nondimensionalizing the equation of 
motion (3.6). By analogy to SII, we define the dimension- 
less coordinate rj = a;/(i>a;)^/^ = x/£{uj) = x {A/ujC,)~^^^ 
and rewrite the solution as 



y(.T,^)-2yo5^K"*/^(^?)} 
and Eq. 1.2 as 

ih{7i) = -a^/i (77) . 

The solutions of (5.7) are of the form 



h{n) 



ce 



(5.6) 



(5.7) 



(5.8) 



where 7 may be any one of the four distinct (complex) 
numbers such that 7^ = — i. These are 

^■=i3e-"/» (j = 1...4) . (5.9) 

The general solution is the sum of these four solutions, 

4 



where 



0.92- 0.38 i . 



(5.10) 



(5.11) 



The unpleasant (but certainly not subtle) remainder of 
the problem is to solve for the four cj 's, given some four 



boundary conditions. At the left end, we enforce not only 
the position at x = but the condition of torquelessness, 
as appropriate for an optical trap, yxx{0) = 0. The right 
end must satisfy the free end boundary conditions, Eq. 
2.7. The Cj derived from these conditions are functions of 
a rescaled polymer length C = L/£{co) and may properly 
be written as Cj (£). 



B. Semiinfinite polymer 

The exact solution for h (77) is presented in Appendix 
C; it simplifies greatly, however, for extreme values of 
L/£{(j). In the limit of an infinitely long polymer, the 
two coefficients Cj for which 7^ has a nonnegative real 
part must be zero, allowing only decaying solutions as 
X ^ 00. 




L/l(«)=l — 



8 I — ' 



_l I I L. 



2 4 6 

r]=x/\{u) 

FIG. 8. Solutions to BHD problem II for filaments of vari- 
ous rescaled lengths C 

The solution consistent with the two left-end boundary 
conditions is 



c-^**" cos {SsV + t^t) 
+ e-^«^ cos (CgT? - iot) 



(5.12) 



where Cg = cos (tt/S) and 5*8 = sin (tt/S). Compare with 
the solution to SII, Eq. 5.2. The semiinfinite solution 
(5.12) is shown at the bottom of Fig. 8 for tot = n27r/6, 
n = 1 . . . 6. In the hydrodynamic case, the solution 5.2 
describes exponentially decaying right-moving traveling 
waves of transverse velocity. In the elastohydrodynamic 
case, the higher-order derivative allows more complicated 
behavior: right and left-moving waves of displacement, 
with different decay rates and velocities. In this case. 
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the right-movers have a slower decay (since 5*8 ~ 0.38 < 
0.92 ~ Cs), and might be expected in some sense to 
dominate over the left- movers. This mechanism will be 
elaborated on in section VD. 



C. Finite polymer 

In the limit of a short or stiff polymer, C <^ 1, we 
may expand the coefficients cj (£) appearing in the exact 
solution to find 

{r]) ~ {zoC/8) [sm{zQr]) - sinh(^o??)] 
+ (1/2) [cos (2:0??) + cosh (-20??)] 
- (3/42;o£) [sin (^0??) + sinh (^0??)] • (5.13) 

The apparently singular result h ^ L^^ is shown illu- 
sory by rewriting r] = aC, a = a;/Le(0,l) and expand- 
ing (5.13) for small £, yielding 



/i« (a) ~ (^1 - 



-a 



+— - (-16 + TOa^ - 70q-' + 21f/ 
1680 



(5.14) 



Equivalently, we may derive this polynomial by truncat- 
ing a series expansion for h'm. a and enforcing the equa- 
tion of motion (5.7) and the boimdary conditions (2.7). 

Using (rj) (in 5.13), the boundary conditions at 
X = and Eq. 5.7 are satisfied exactly, whereas the 
boundary conditions at x = L are satisfied to order 
O ('C^). Using (a) (in 5.14), all four boundary con- 
ditions are satisfied exactly, whereas Eq. 5.7 is solved to 
order O (£4). 

The exact solution is shown in Fig. 8 for C ~ 1,2,4, 
and 00 and Lot = n27r/6, n — 1 ... 6. Note the existence 
of a pivot point at x = 2L/3 as £ — * O. This behavior 
is described by the 0{£^) term in Eq. 5.14: as £ ^ 0, 
the polymer acts as a rigid rod. As a consequence, it 
is impossible to tell if a movie of such a polymer is be- 
ing played forward or backward. Indeed, this is a fila- 
mentous version of the famous "one-armed swimmer" or 
"scallop" example illustrating the lack of net propulsion 
for rigid objects executing time-reversible motions in low 
Reynolds number flow [4, 28]. 



D. Propulsive force 

Problem II and its associated experiment are suffi- 
ciently reminiscent of flagellar hydrodynamics to moti- 
vate a calculation of the propulsive force F generated in 
the X direction by the wiggling. This can be done by inte- 
grating ff , the force exerted by the polymer on the fluid, 
along the length of the filament. We then contract this 
instantaneous total force with ex and average over one 
period. This force is equal and opposite the propulsive 
force exerted by the fluid on the polymer. 



Noting that the force per unit length in Eq. 
total derivative. 



ff = Ads { '«sn — TT't^t 



2.5 is a 



(5.15) 



and recalling the boundary conditions imposed on k and 
Kg, we have 



/ 



ds fs • Bx = — ^Bx • yKsii — (5.16) 

= Aex • (K«n) (s = 0) (5.17) 
= Aks sin 9{s = 0) . 



This is geometrically exa(;t. We now wish to calculate 
the time average over one period, which we denote by 
(•••). Within the linearized solution, Kg sin ^ ~ UxxxUx- 
Recalling the expression for y in Eq. 5.6, and expressing 
the answer in terms of the rescaled variable rj, we obtain 



F = A{yxyxxx) 
Ayl 



VoO^r 

4\/2 ^ V^(^) 



(5.18) 
(5.19) 

(5.20) 



where * indicates complex conjugation, £{u) is the char- 
acteristic length, and T is a scaling function conveniently 
normalized (see below). 

The exact solution to EHDII given in Appendix C can 
be used to calculate the function T for all values of the 
polymer length. The asymptotic behavior as £ —> 00 is 



T(£) ^ 1 + 4e-2^«^ sin {2CsC) . 



(5.21) 



When the length is short compared to the characteristic 
length, so the polymer flexes very little. 



3360 



3360 A 



(5.22) 



As Fig. 9 illustrates, the short-length approximation 
(5.22) shows good agreement with the exact solution for 
£ < 3, as does the large-£ approximation (5.21) for 
C ^ 3. The approach to the asymptotic limit is oscil- 
latory, with a maximum near £ ~ 4, the value at which 
acquires its flrst root, and a local minimum near 
£ ^ 6, the the value at which acquires its sec- 

ond root. The unexpected local maximum indicates that 
there is an optimal combination of A, w, and a finite L. 
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^"fcio '^*y^^^=''^* dxyix) (5.24) 

or, noting that / dx y{x, t) is simply the area A{t) under 
the curve y{x,t), and that the slope at the left is to first 
order simply the tangent angle 



1^ 




FIG. 9. Scaling function T for propulsive force. The large 
jC expansion is plotted for £ > 2, and the small-£ solution is 
plotted for £ < 3.5. 

Inserting typical numbers ft-om the experiment (in cgs) , 



10-2 



2/0 ~ 4 X 10"* 

277, # 



L ^ ' " = 10^ 

<^ 7 X 10-7 



(5.23) 



we find F{oo) ~ 2 x lO-'' dynes = 3 x 10-^ pN. For a 
trap stiffness of ^ .02 pN/nm, this would induce a dis- 
placement of 1.5 nm, at the lower limit of experimental 
observation. The production and measurement of propul- 
sive force by an artificial fiagellum was attempted by G. 
I. Taylor [32] using a glycerine-filled tub to mimic the 
low Reynolds numbers found in vivo. Taylor struggled 
to drive the fiagellum without inducing unwanted torque 
or disturbing the flow, a difiiculty obviated by the use of 
optical traps. 

Returning to the asymptotic expressions for h derived 
in the sections VB and VC, we observe a pleasant ac- 
cordance with the qualitative features of Fig. 9. In the 
semiinfinite case, we noted the presence of right- and left- 
movers, with right-moving waves of displacement exhibit- 
ing slower decay. Such a dominance accounts for the 
nonzero propulsive force in the £ — » oo limit, where a 
net propulsion to the left survives. In the £ ^ case, we 
recovered a shape which asymptotes to a pivoting rigid 
rod, not unlike a one-armed swimmer. As we expect 
from life at low Reynolds number [4, 28], such a motion, 
invariant under t —t, can produce no net propulsion. 

As further illustration of the relationship between 
low- Reynolds- number swimming and cyclic motions, we 
observe that the lowest-order expression for the time- 
averaged force is equal to 



F = 



27r 



dt 6*0 



dA 

dt 



9odA. 



(5.25) 
(5.26) 



This result can be interpreted quite simply: the propul- 
sive force results from pushing aside some volume (or 
in two dimensions an area) of fluid, projected in the di- 
rection of propulsion ex an amount proportional to 9o- 
Note that had we been interested in the propulsion in 
the transverse (ey) direction, the Oq would not appear, 
leaving the absence of net forcing: F (x f dA — 0, as we 
would expect. 

The net force, then, is the area enclosed by a trajectory 
in A— 9o space during some cyclic motion. This represen- 
tation is independent of the particular motion exhibited, 
although we have here considered simple periodic motion, 
for which the trajectory is always an ellipse. As £ — > 0, 
the elliptical trajectory thins to a straight line, encloses 
no area, and thus produces no force, as commented on 
above. 

This representation makes clear that in an inertialess 
world, net motion is principally geometric in origin rather 
than dynamic. [36] In a manner analogous to the impor- 
tance of path rather than kinetics in generating net work 
in a Carnot diagram, we see that we can remove time en- 
tirely from the expression and consider instead a path in 
a low-dimensional projection of the infinite-dimensional 
shape space. 



VI. EXPERIMENTS ON ACTIN 

As mentioned in section V, the EHDII experiment was 
performed and the data compared to the solution of Eq. 
5.7. In this way we were able to confirm the results of 
the analysis and investigate a new method for measuring 
biopolymer bending moduli. 



A. Materials and methods 

Actin from white leghorn chicken breast was pu- 
rified after a published procedure [27]. Filamentous 
actin (F-actin) was fiuorescently labeled with rhodamine- 
phalloidin (R-415, Molecular Probes, Inc.). Anti-actin 

antibodies (Sigma) were covalently coupled through car- 
bodiimide to fiuorescent polystyrene beads (CX, Duke 
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Scientific Corporation), of diameter 1 mm, after the rec- 
ommended protocol of Duke. Experiments were per- 
formed in actin suspension buffer containing 25 mM imi- 
dazole, 25 mM KCl, 5 mM b-mercaptoethanol at pH 7.65. 
To avoid photobleaching during observation, Img/ml glu- 
cose, 33 units/ml glucose oxidase, 50 units/ml catalase 
were added to the suspension buffer. 

In order to prevent actin filaments from sticking to 
the glass surfaces, slides were coated with bovine serum 
albumin (BSA). Observation chambers were sealed with 
nail polish and filled by capillary action. Beads and F- 
actin coupled within the cell; the ratio beads:F-actin was 
adjusted to have one or two beads per filament. Beads 
bound to various locations along filaments. Among these 
filaments, we chose bead-ended actin filaments for the 
experiment. 

Observations were made on a Zeiss microscope (Ax- 
iovert 135) equipped with a 100 W mercury lamp and a 
standard filter set (XF37, Omega optical). To prolong 
observation time, the excitation light intensity was re- 
duced by inserting neutral density filters. Fluorescent 
images were taken through a 4X TV tube, via an im- 
age intensifier (Hamamatsu) followed by a CCD camera 
(XC-77, Sony). Images were recorded with an S-VHS 
recorder. Images corresponded to an overall screen size 
of 44 mm in horizontal dimension leading to a value of 
0.06mm/pixel after digitization. Fig. 7 shows the ex- 
perimental setup, and Fig. 10 shows a series of typical 
images. The smooth curves indicate fit solutions to Eq. 
5.7. 




2 4 6 8 
X [/i.m] 



FIG. 10. Series of typical images of driven actin filaments 
in EHDII. Also shown are fits to the solution of Eq. 5.7. 

The optical tweezer was made by focusing a 0.5 W 
Nd:YAG (Spectra) laser beam through a Zeiss 63X, 1.4 
numerical aperture Plan-Apochromat microscope objec- 
tive. A mirror mounted on a galvanometer controlled by 



a function generator served to move trapped beads in the 
focal plane (Fig. 7). We imposed sinusoidal driving, with 
frequencies and amplitudes ranging, respectively, from 
0.1 Hz to 6 Hz and from 5 jim to 10 /um. 

Up to 10 images per oscillation period were digitized 
with a Scion frame grabber card and a Power Macin- 
tosh; they were subsequently analyzed using NIH-Image 
software. Each image was retraced by mouse and its co- 
ordinates determined. Only those sequences for which 
filaments remained in the focal plane for an oscillation 
period were retained for data analysis. 

B. Experimental results 

Knowing the amplitude of the driving of the bead (yo) 
and the frequency (w), and reading off the projected 
length (L) and the phase (cut) directly from the images, 
we were left with a one-parameter fit of the images to 
the solution of Eq. 5.7, varying only £{iL>) to minimize 
X^. We can then observe the dependence of £ on w, as 
illustrated in Fig. 11. The variation in error bars can be 
attributed to the widely-varying numer of images taken 
at different frequencies. 

Comparing with the earlier analysis (c/. section V A), 
we can extract from this scaling a measurement of the 
persistence length. Fitting to 

^(o;) = (M:y''ii/4^-V4, (6.1) 
we determine Lp to be 7.37 ± 0.25/im. 



2 



1 






1 


- ] 


\^ 






L 
1 


p=7.37±0.25/xm 


1 , . . , 1 . , . 





10 20 30 40 

<a [Hz] 



13 



FIG. 11. The characteristic length scale £{u)) vs. frequency 
uj. The smooth curve is a fit to Eq. 6.1 for Lp = T.ST/um. The 
scaling of the characteristic length with the fourth root of time 
suggests that actin is well-described by Eq. 1.2 and thus is a 
semiflexible polymer, with a scale-independent elasticity. 

C. Discussion 

There are a few limitations with this realization of 
the experiment which, upon correction, will improve this 
technique and make the data more conclusive. An ob- 
vious mechanism for improving the error bars is to ac- 
cumulate more data. Image-taking was entirely manual, 
with hundreds of student-hours spent looking for usable 
images. An automation of this process would clearly be 
advantageous and improve the low statistics used here. 

With careful control of the timing, images of equal 
phase can be superposed to average out thermal fluctua- 
tions or experimental variation in the images before fits 
are performed. 

Most importantly, since our aim was to verify the plau- 
sibility of the experiment, we did not limit ourselves to 
small-amplitude wiggling, thus leaving the realm of va- 
lidity of the small-j/x approximation. We can estimate 
the error due to such driving by looking at the relevant 
terms from the geometrically exact equation of motion: 

yt ~ I 1 (R o\ 

If we wish to approximate this with Eq. 3.1, we are 
measuring an "effective" or "i/" , where 

V ~ (6.3) 

where the brackets indicate averaging over the data, and 
thus the true A will be underestimated by a factor of 
{{1 +y^)'^), which is always greater than unity. Inspect- 
ing Fig. 10, we see that there are data for which is 
not necessarily small. For this reason, the data we col- 
lected can only put a lower bound on A. Wc anticipate 
that the true value may be greater by a factor up to 
1.5; future experiments clearly should employ smaller- 
amplitude driving. 



experiment. Similarly, in an attempt to make the anal- 
ysis as clear as possible, wc have omitted from Eq. 3.1 
the background flow due to the trapped bead. 

Some amount of discussion has been entertained in 
the biophysical community about the possible scale- or 
time-dependence of the elasticity of biopolymers, includ- 
ing both actin [16] and microtubules [18]. One of the 
more powerful features of this technique is that, since it 
involves a controlled dynamic, specific scales and frequen- 
cies can be investigated to attempt a spectroscopy of elas- 
ticity. Equation 3.1 can be extended without difficulty to 
include a characteristic relaxation time t or a continuum 
of times, in an attempt to model a characteristic rate of 
bond-breaking in the presence of bending. Similarly, one 
can include additional bending moduli which depend on 
higher-order derivatives. For oscillatory motion, occur- 
rences of dt are simply replaced by dt + 1/t or iw + 1/r. 
Including higher-order derivative terms simply results in 
replacing Ay2x with Ay2x + By4x in the bending energy 
and Ay^x with Ay/^x + By^^ in the equations of motion. 
We then may recover such an equation as 

C{dt + -)y = -Ayxxxx - By^x, (6.4) 

T 

which may be solved in the manner of Eq. 1.2. This 
more general expression makes possible the experimental 
confirmation or refutation of such hypothesized mecha- 
nisms. 



VII. COMMENTS ON PREVIOUS 
EXPERIMENTS 

A recent pair of elastohydrodynamic experiments in- 
volving microtubules [10], brought to our attention as 
the original version of this paper was being completed, 
provides an excellent opportunity to apply the spirit of 
analysis which we have developed for EHDI. In both, the 
crucial experimental observable is the motion of the free 
end of a microtubule; the analysis must then relate this 
motion to the bending modulus A. We will first investi- 
gate the analysis appropriate to a special case of EHDI 
in which a polymer relaxes to a straight configuration in 
the absence of a driving fiow. We then investigate a more 
complicated experiment in which an optical trap induces 
a force in the middle of the polymer. 



D. Suggestions for further analysis 

We also anticipate that a more accurate treatment of 
the geometry and hydrodynamics would refine the tech- 
nique. The true geometry is nonlinear and hydrodynam- 
ics nonlocal, but neither intractable and both amenable 
to numerics. The geometrically exact, intrinsic formu- 
lation involves some enjoyable mathematics of curve dy- 
namics, [3] whereas the linearized treatment presented 
herein is more illustrative and more easily connected with 



A. The Simple EHDI experiment 

In the first experiment [10], a microtubule is initially 
clamped at the left to an axoneme and trapped directly 
at the free end; the polymer is bent out of its mechanical 
equilibrium configuration. When the trap is shut off, 
the polymer relaxes back to the straight shape in a way 
which we may describe as before: the initial condition 
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is projected onto the appropriate space, in which each 
mode decays independently. 




(b) 

FIG. 12. EHDI experiments, (a) a simple special case (b) 
end-driving as well as driving via a point force. 

What remains then is merely to determine the initial 
data: the shape of a biopolymer clamped at one end and 
held by a trap at the other. We may derive this from 
the (geometrically exact) general equations of force and 
torque balance for the elastica: 

N« + fe = (7.1) 
m, + (r X N)^ + r X fe = (7.2) 

Here N is a contact force (not the unit normal), m is a 
bending moment, and fe is some external force. While 
these equations may appear unfamiliar, they in fact have 
a long history as the classical forms of the equations of 
equilibrium for the special theory of Cosserat rods [2]. 
Combining the equations into one, 

nis = t X / dafe (cr) , (7.3) 
Jo 

where t = is the unit tangent. In two dimensions, 
cross products are scalars, and mg becomes the scalar 

We now include the defining constitutive relation of 



linear elasticity:^ m = AOs- Furthermore, in considering 
point external forces (the axoneme and the trap) that act 
at the left and right ends, we write 



fe{s)=Fo6{s) + Fi6{s-L) . (7.4) 

We assume that there is no net force, so Fq = — Fi. 

Noting that t = (cos ^, sin ^) we rewrite Eq. 7.3 (for 
< s < i) as 

Ae^s =- Fy cos e + sin 6 (7.5) 



where F^^ yy are the components of the force exerted by 
the trap, rather than by the axoneme. We expect Fx < 
and Fy > 0. Eq. 7.4 is geometrically exact and can be 
solved in terms of elliptic integrals. 

Note that there is no reason to assume the trap ex- 
erts a force only in the By direction - the oft-quoted 
"cantilcvered beam" problem from introductory civil en- 
gineering texts. We show below, however, that the ex 
component is higher order in d and will be ignored in the 
linearized treatment presented. Linearizing Eq. 7.5 for 
small 9 obtains 

Ae,, = -Fy + 9Fx, (7.6) 

the windy pendTilum equation [35]. Assuming the poly- 
mer is homogeneous along its arclength, we note that 
there is no energy cost upon moving the trap along the 
ax:is of the polymer. Therefore the trap can exert no force 
in the tangential direction. We may write this constraint 

= F -ic^ Fx + Fy^ . (7.7) 

Accordingly, F^/Fy ~ —dy/dx, and Eq. 7.6 becomes 

Ayx.x^-Fy{l + 0{yl)) . (7.8) 

We now must consider the boundary conditions. The 
axoneme clamps the left end of the polymer, thus 
y (0) = yx (0) = 0. Since there is no energy cost to ro- 
tating a polymer held in an optical trap, there should be 
no bending moment, implying y^x (L) = 0. Given these 
three boundary conditions, the solution to (7.8) is 

y{a:) = l{3a'-a') , y = y(L) = ^ , (7.9) 

where again a = x/L. 

Now that we have our initial data, we project it onto 
the appropriate function space in which the dynamics are 
simply exponential relaxation. In keeping with the exper- 
iment we seek here to model, we focus on the motion of 



^Hcrc, "linear" refers to the curvature in the constitutive 
relation, not necessarily in the equations of motion, in which 
the geometry is often the source of nonlinearity. 
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the free end, whose rate of relaxation provides a direct 
means of measuring the bending moduhis A. Given the 
dramatic increase in relaxation rates for each subsequent 
mode (c/. Eq. 4.14), we expect only the lowest mode 
to be relevant beyond negligible initial times. Moreover, 
inspecting Fig. 6, we observe that the eigenfunction Wki 
well-approximates the normalized third-order polynomial 
which describes the initial data, whereas higher modes 
more closely resemble Fourier modes. This close agree- 
ment, rather surprising from a sum of trigonometric and 
hyperbolic trigonometric functions, leads to the domi- 
nance of the projection of initial data onto the first mode; 
specifically 

Vk, (0) ~ 0.971 y; yu, (0) ~ -0.0247 y , (7.10) 

where we define yfe„ = /g^ da yWk^ , in analogy to Eq. 
4.14. Inserting typical numbers (A 5 x 10~^^ dynes 
cm^, L ~ lOx lO^m, and the drag coefficient from Eq. 3.3 
with d ~ 20nm and /i ^ .01 erg s/cm"*), we see that for 
times beyond 0.01s, the amplitudes of subsequent modes 
are at most 1% that of the first mode. The shape is 
then described by y (x) ~ Wfei {x) yui (t) and decays as 
e~'"*, with r = Ak'l/C,L'^. In the model accompanying the 
experiment, it was assumed that the shape was described 
by a single decaying mode for all times up to t = 0. 
As we have shown, Wki fortunately well-approximates 
the initial data such that this introduces only an error 
0(10-2). 

Within these approximations, the free end decays as 



y{L,t) 
y{L,o) 



0.971Wfc, (L)e 



-rt 



-rt 



(7.11) 



so that a measurement of the e-folding time yields r and 
hence A. It is now clear that the identity of fci, the 
solution to the transcendental equation (4.17), is crucial, 
as it is physically manifested in the decay rate. 

The model in [10] is based first on the computation of 
the deflection yi of the free end of the clamped elastica 
experiencing a flow linearly increasing in x and constant 
in time, with a drag coefficient C. Such a fiow would be 
appropriate for a rigidly rotating rod with constant an- 
gular velocity, rather than for a bent filament. Second, to 
find the decay rate in [10], it is noted that the exponen- 
tial relaxation of the tip must correspond to some first- 
order differential equation. Although this equation is not 
stated, we must assume it to be yt{t) = ■~{v/yi)y{t), 
(where v is the maximum flow velocity used in the com- 
putation of yi) in order to recover the reported decay 
rate. 



r = 



120 "yl" 



(7.12) 



where "A" refers to the value of A which would be ex- 
tracted from data, using the simplifications described 
above. 



A careful analysis of the drag coefficient in slender- 
body hydrodynamics employs a matched asymptotics for 
the fiuid velocity and is dependent on the shape of the 
slender body. Treating a microtubule as a cylinder, the 
appropriate drag coefficient is [6] 



c = 



47r/z 



ln(L/rf) + 21n2- 



(7.13) 



Unfortunately, the drag coefficient used in [10] is that 
appropriate to tangential rather than normal flow, with 
numerator 27r/U (quoting an earlier misquotation of Doi 
and Edwards [9, 33]) and further suffers from the re- 
placement of the constant terms in the denominator of 
Eq. 7.13 with - In 2. 

We may now compare the results of a differential 
equation- motivated analysis with the model of [10]. 
Given some measured e-folding time U, the eigenmode 
analysis yields the bending modulus A as 



(7.14) 



whereas the rodlike treatment implies (from Eq. 7.12) 



C,L^ 11 
IT 120' 



(7.15) 



Inserting typical numbers from the experiment, L 
10/xm, d ^ .02/^m, we see that 



"A" 



1.13 



0.709, 



(7.16) 



a systematic underestimate beyond the uncertainties of 
experiments. 



B. The discontinuous EHDI experiment 

A more complicated example of EHDI which was also 
conducted includes driving the axoneme, so that we must 
incorporate an inhomogeneous boundary condition, and 
point forcing by an optical trap [10]. In this experiment, 
the axoneme (attached to the cover slip) is moved with 
constant velocity between two extreme positions. During 
this motion, the polymer is constrained by the optical 
trap to pass through some intermediate point {xp,yp). 
From the position and velocity of the free end it is pos- 
sible to determine the bending modulus. Enumerating 
all the relevant forces, we consider the force at the left 
end due to the axoneme, F^, the point force due to the 
trap Fp, and the force per unit length due to drag, —C,ht 
{cf. Fig. 12). We insert these terms into the equations 
of force and torque balance (7.1,7.2) to find the equations 
of motion. 
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1. Explicit declaration of forces 

The sum of the external forces per unit length can be 
expressed as 



fe -.= ¥^5 {s - Sp) + ¥J {s) + id 



(7.17) 



where is the drag (Eq. 3.1). 

Integrating over the delta functions, and taking the 
axoneme to be located at s = 0, the geometrically ex- 
act governing equation (7.3) becomes the fiercely com- 
plicated integro- differential equation in 6{s, t) 



AOs 



t X 



(7.18) 



+C, da (fin -|- (3ii) (^J^ da'n9t - u 



which greatly simplifies for small ■ Here u is the back- 
ground flow velocity, in the frame of the trap, due to 
the motion of the coverslip. This is a constant for the 
problem, which we define to be —VcSy so that Vc > 0. 
Note that the axoneme, attached to the coverslip, also 
moves with velocity —VcGy. We see that point forcing in 
the middle of a polymer introduces a qualitative change 
in the dynamic: a Heaviside function in the equation of 
motion. 



2. Linearized geometry 

We adopt an expansion in 9, keeping only the first- 
order terms. The drag term, as in Eq. 3.5, simplifies to 
Cey{yt+Vc). 

Noting, as in Eq. 7.7, the constraint that the force on 
a polymer due to an optical trap must have no tangential 
component, the y-componcnt of the linearized equation 
reduces to A0ss — Fa — FpQ {s — Sp) — C(/'' yt+Vc), which 
upon differentiation and linearization implies 



CiVt + Vc) = -Ay^ 



FpS {x - Xp) . 



(7.19) 



This is our working equation, obtained from the lin- 
earized second derivative of the equation of net torque- 

lessncss. 

Since we do not know the magnitude of Fp a priori, we 
must perform a matching of (x) , the curve described 
by the anchored end, and y^ (x), the curve described by 
the free end, which solve 



.^y{a,/} 

yxxxx 



(7.20) 



subject to matching conditions at the point of forcing. 
Horrifying though this may sound, it is an excellent op- 
portunity to apply the simple ideas developed in earlier 
sections. It is also, of course, a chance to treat in a thor- 
ough way the analysis of a biopolymer subject both to 
drag and to micromanipulation via some point force and 



boundary condition. Moreover, if we wish to claim that 
the PDE description is the appropriate analysis for single 
biopolymer dynamics, we must be able to apply it even 
to such an awkward case. We will separate the solutions 
into the homogeneous and the particular, and in a pro- 
cedure that is now familiar, construct the appropriate 
function space in which the dynamic is simple. 



3. Matching and boundary conditions 

Now that wc have two fourth-order equations of mo- 
tion, we must specify eight matching and boundary con- 
ditions. Inspecting the equation of motion (Eq. 7.19), we 
see it supports a discontinuity in yxxx', however, y^x, de- 
scribed by an integral over a Heaviside function, is contin- 
uous, as are y and yx- Moreover, if we wish to describe an 
experiment in which the filament position is constrained 
at the point of forcing Xp, we have not only the matching 
condition y°'{xp) = y^{xp), but the stronger condition 
y°'{xp) = y^{xp) = yp-, here we choose yp to be with- 
out loss of generality. We are thus describing a polymer 
pinned at a certain intermediate point along the curve 
by an optical trap, while the right and left sides perform 
some coupled motion. 

At the free end, wc impose forcclcssncss and torque- 
lessness, stated linearly as yxxx {L) = and yxx (L) = 0, 
respectively. At the anchored end, we describe a poly- 
mer clamped to the axoneme position yA(t), thus y (0) = 
yA{t); Z/x (0) = 0. The above-listed boundary condi- 
tions constrain 4 of the arbitrary constants, the remain- 
ing matching equations constitute the remaining 4, and 
we may thus completely specify the solution. 



4- Construction of the inhomogeneous solution 

The governing equation is linear, allowing us to sep- 
arate y into two separate solutions of Eq. 7.20: y = 

yo{h + g), where yo is some typical length scale. Wc 
choose g {h) to satisfy the (in)homogeneous boundary 
conditions and equation of motion. Compare this with 
the first example solution of EHDI (section IV B), in 
which h was chosen to satisfy an inhomogeneous equa- 
tion of motion but homogeneous boundary conditions. 

As described above, the axoneme moves with con- 
stant velocity Vc from yA = A to yA = —A: yA (t) = 
—Vet, {—A/vc < t < A/i'c). Wc must now merely solve 
for the inhomogeneous solution, expressed as 



hi-J} (a;,t) = ^c{«'/> {t)x^. 



(7.21) 



n=0 



However, as in the derivation of Eq. 4.24, wc expect this 
solution to depend only linearly on the driving and thus 
linearly on time. Constraining dfcl^''^^ = 0, enforcing 
the 8 matching and boundary conditions, and respecting 
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the relationship between 9tc„ and c„+4 dictated by Eq. 
7.20, we completely specify the solution. 



5. Polynomial solution 

We here quote the solutions for In order to 

make explicit the qualitative behavior in limiting cases, 
and to write the solutions as compactly as possible, we 
employ three nondimensionalized variables: 



a = 



X 



a 



■ = t 



(7.22) 



The coordinates {a, a} G (0, 1) measure distance on the 
left from the axoneme and on the right from the trap, 

while T is the time rescaled by the characteristic elasto- 
hydrodynamic time for the anchored section. The fact 

appears explicitly in 



that Xp, rather than I 



our choice of definition of r is reflected in the equation of 
motion, in that the functions now solve slightly diflferent 
equations for the two sides: 



yo A 
yo 



(7.23) 



The natural choice for yo is clearly yo = VcCXp/A, where- 
upon we are left with the dimensionless equations of mo- 
tion 



K 

hl + l 



1 = 



where we have introduced the ratio of lengths 



P = 



L - .-r„ 



I 



(7.24) 



(7.25) 



which describes the location of the point of forcing Xp 
relative to the "extra" length I. AsXp nears the anchored 
end X = 0, p ^ 0, and as Xp nears the free end x = L, 
Xp/ {L — Xp) = p ^ oo. 
The expressions are then: 



1 



80 



1 4 

— pa 

24' 



6 8/ 



ZP 



1 3 ^ 
^P+^Pr]a 



(7.26) 



6. Decaying modes 

As described before, the solution h to the inhomoge- 
neous boundary conditions will describe the long-time be- 
havior remaining as transients decay exponentially. We 
now turn our attention to the transient g, which must 
satisfy the following homogeneous boundary and match- 
ing conditions: 



5(0) = 0, .9.(0) = 0, 
gxx{L) = 0, gxxx{L) = Q 
g (xp) = 0, [go: (xp)] = 0, [fif^^x (xp)] = 0, 



(7.28) 



where the brackets indicate discontinuity. 

The fact that all the boundary conditions are 0- valued 
means we stand a chance of constructing a self-adjoint 
operator, consistent with these conditions, from the rele- 
vant differential operator: d^. To do so, we left-multiply 
the equation of motion by an as-yet arbitrary function 
Wk (x) and integrate over the entire domain. Since h 
is constructed to solve the (linear) equation of motion 
(7.19), g must as well, and we derive an equation of mo- 
tion** for the quantities gk = Jq dxWk{x)g{x): 



Irdtgk 



Jo 



dxWk {x) gx 



(7.29) 



We then wish to integrate the RHS by parts. However, 
we must admit the possibility that g supports a disconti- 
nuity in its higher-order derivatives, i.e., [gxxx (xp)] ^ 
0. For this reason, we define {5" (x, t) , (x)} and 

{g^ {x,t) ,yvl {x)}, as in Eq. 7.20, defined on the an- 
chored and free sections, respectively. If we choose W 
to obey the same homogeneous boundary conditions as 
g, all surface terms from the integration by parts in Eq. 
7.29 vanish except 



- a,W,"a^5« {xp) + dxWidlgf {xp) 
+dlWtdxg'' {xp) - dlwidxgf {xp) 



(7.30) 



which vanish upon choosing W to obey the same match- 
ing conditions as g. The remaining equation is 



- \dtgk 

V 



fdtWi 







(7.31) 



,{«■/} 



to obey the eigenvalue condition 

Since we have constructed a self- 



We now choose Wl 

adjoint operator, we can be sure the eigenvalues are real 
and positive (c/. Appendix B). We are left with 



llo 1 1 3q\9 

p + -p 4 h -p T 

560^ 8^^ 8 2^^ 



(7.27) 



**The discontinuity in yxxx is shared between h and g, so we 
might for the sake of explicitness include a term on the RHS 
of (7.29) proportional to Wkixp). However, wc will choose 
Wfc(a;p) = and thus this will not complicate the behavior of 
9k {t). 



18 



pXp pL 
Jo Jxp 

= vk^gk, (7.33) 

the solution of which is, as before, gk {t) = gu (0) q-"^"^^ 
for A; 7^ 0. Fortunately, the boundary and matching con- 
ditions do not admit a solution to gxxxx = 0, and we need 
not consider a mode. 

Given some initial condition y{x, 0), we project it onto 
this strange cigcnspacc spanned by {Wfe}. Wc then con- 
struct g {x, t) for all later times in a standard Green's 
function way, 

g{x,t)= f'^g{x,x';t){y^^-h{x',Q)) (7.34) 
where 

g = J2 e-''='*Wfc (x) Wk ix') , (7.35) 

k 

and we see that all modes die. 
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FIG. 13. p — q plane showing roots of the solvability eon- 
dition. Here p = kxp and q = kl = k{L — Xp). The geometry 
of the experiment dictates p, as described in the text. Verti- 
cal and horizontal lines correspond to solutions of F{p) = 
and H{q) — 0, respectively. The diagonal line indicates the 
solution to Xp/L = .8. 

7. Construction ofWk; solvability condition 

We now solve for the coimtably infinite sets {k} and 
{Wfc}. The general solution of d^wl"'^'"^ = k'^wl"'^^ is 

y^{a,/} ^ ^{aj} ^.^ ^^^^ ^ ^{aj} ^^^^ 

+4"'^'^ sinh {kx) + ai""'^^ cosh (kx) , (7.36) 

where the 8 arbitrary constants will solve the 4 boundary 
and 4 matching conditions. The insertion and elimina- 
tion of these constants is not a joyful task and will be 



omitted here. The most important fact is that the set of 
8 equations for 8 unknowns can be written as a matrix 
of trigonometric functions (as in Eq. 4.16) annihilating 
the vector a of unknowns. The 0-valuedness of the con- 
ditions dictates that the determinant V of this matrix be 

0. This solvability condition is written explicitly in ap- 
pendix D; the solutions, graphically constructed in Fig. 
13, determine the allowed values of k given some fixed 
ratio p. 

Curiously, the complicated solvability condition can be 
expressed compactly as the separable equation 

V = dpF ip) H{q) + F ip) dgH (q) = 0, (7.37) 

where F{p) = cos(p) cosh(j3) — 1, H{q) = cos(g) cosh(g) -t- 

1, p = kXp, and q = kl. This differential relation de- 
scribes the motion along each of the branches shown in 
the figure as p varies, each branch indexed by arbitrary 
constants introduced upon integrations of Eq. 7.37, and 
separated by the singularity lines {H = 0, = 0} at 
which the differential equation is not invertible. The ge- 
ometry chosen by the experimenter dictates Xp,l, and 
therefore Xp/l = p/q = p. Inspecting the figure, we see 
that we choose a set of modes by drawing a line of slope 
p through the origin; each intersection with a branch cor- 
responds to one mode. 

A further curiosity is that the each of the four equa- 
tions {F, H, Fp, Hq} = is itself a separate solvability 
condition associated with a separate experimental geom- 
etry. In the language of Appendix D, we may rewrite the 
condition V = 0, with the first letter below indicating 
the boundary condition for the left end of one side, and 
the second for the condition at the right, 

V = {cf){hf) + {ch){ff). (7.38) 




0.0 0.2 0.4 0.6 O.B 1.0 

x/L 

FIG. 14. The first three modes for the discontinuous 
EHDI experiment. The circle indicates the pinning point 
Xp/L = p/{p + q). 
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Otherwise stated, the total solvability condition is an 
average of those for (i) the left side clamped at and free 
at Xp, with the right side hinged at Xp and free at L, and 
(ii) the left side clamped at and hinged at Xp, with the 
right side free at Xp and free at L. 



8. Initial data 

Armed with a set of decaying modes, we may solve the 
equation for all times. We now attempt, as before, to 
construct a polynomial solution which describes the ini- 
tial data. The polymer sits at rest with y(0) = A and is 
subject to the stated matching and boundary conditions. 
Noting that the shape must be described by polynomi- 
als of less than fourth order, since we wish to describe 
an elastica experiencing no forces, we find for the new 
polynomial y 



r («) = A ( 1 - + 



{a) = -A 



So- 



(7.39) 
(7.40) 



The final configuration, after the axoneme comes to rest 
and all transients have died, will be —y: a third-order 
polynomial on the left and a straight line on the right. A 
pleasant fact is that this polynomial can also be derived 
by taking the limit as t;c — * of yoh^"''^^. 



9. Comparison with experiment 

We now wish to use this information to arrive at a mea- 
surement of A. Hoping to verify the plausibility of this 
analysis, we compare with the results and accompanying 
model published with the experiment. 

The model presented in [10] can be summarized as fol- 
lows. First, the forces are calculated on a slender body 
subject to (i) a constant flow and (ii) a flow linearly in- 
creasing from at the origin to some Vm at the end. The 
linearly-increasing flow describes that experienced by a 
rigidly rotating rod, in hopes that the force experienced 
by the actual (curved) filament, itself neither a straight 
rod, nor moving with constant velocity in time, is well- 
approximated. Curiously, different drag coefficients are 
used for these two forces. 

The forc:es are then \iscd to compiite displacements, 
using the results for the elastica clamped at the origin, a 
condition which is unfortunately consistent with neither 
the assumption of a rodlike shape nor with the experi- 
mental geometry, in which the free side of the elastica is 
neither clamped (j/j, = 0, as was assumed for the calcu- 
lation of displacement) nor hinged {yxx = 0, as might be 
considered appropriate for a straight rod) ; it is in fact ex- 
periencing a torque at Xp due to the side of the filament 
between the trap and the axoneme. 



Summing (i) the displacement yc{vc]Ci) of a clamped 
elastica in constant flow (using the first drag coefficient) 
and (ii) the displacement yi{vm', C2) of an elastica subject 
to the linearly increasing flow it would have experienced 
were it a rigid rod (using the second drag coefficient), 
the total deflection ys from a straight line is obtained 
(see Eq. 20 of [10]): 



ys = 



■7rri{llvm + 30vc)l'^ 
60Aln{2l/d) 



(7.41) 



It is then assumed that this displacement equals the de- 
flection which would have been experienced had the fila- 
ment not been initially horizontal, but rather constrained 
to some nonzero slope at the trap, specifically, that of the 
t ^ 00 solution for the shape, —y. 

This deflection yg, from the right tip position at t = 00 
to that a.t t = A/vc, when the axoneme halts, is the first 
experimental observable. The second is the velocity at 
this tip at t = A/uc, which is then equated with the 
maximum velocity in the expression for yi. Dividing 
ys by the weighted sum of velocities appearing in the nu- 
merator, one obtains the combination of two observables 
{ys,Vm) and one experimental parameter {vc) which, ac- 
cording to this model, is equivalent to a simple quotient 
with units of time and in which A appears explicitly. 



T. 



ys 



iittt/;^ 



(llw„ + 30vc) M"60 ln{2l/d) 



(7.42) 



where "A" indicates the rigidity which would be ex- 
tracted from the data using this model. Note that the 
expression is a function only of Z, the length of the free 
segment of the polymer. In this model, I is taken to be 
"the hydro dynamically relevant length" [10] . 

Returning to the PDE-trcatmcnt of the problem, we 
see that the solutions to Eq. 7.20 for the long-time poly- 
nomial shape relate the velocity of the free end to the 
constant stage velocity as 



3 / 



2 Xj) 



(7.43) 



This simplification allows us to solve for the ratio in 
(7.41) of observables and parameters in terms of the poly- 
nomials in Eqs. 7.40 and 7.26. 



T 



(1^ 



- 16 -h 210 p2 + 420 p -I- 231) 



2520y4 



(ll + 20p) 



(7.44) 



where A indicates the value of the bending modulus 

which one would calculate had one used the same data 
and thus the same quotient T. 

The ratio A/ "A" is a simple expression and allows us to 
compare by what factor this differential equation-based 
analysis differs from the published model given some ex- 
perimental data T{1). Equating the two expressions for 
T and isolating A/ "A" , we find 
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A 

"]4" 



2(1+P) 



210p2 - iGp-^ 
231 + 420p 



(7.45) 



where we have moved the 0(1) constants into the argu- 
ment of the logarithm for compactness. The dominant 
behavior is captured by the polynomial dependence on p, 
but we can see that the function will be ^ 2, indicating a 
systematic underestimate by the published model of the 
bending modulus. A plot of the ratio, with typical values 
taken from experiment, is shown in Fig. 15. 
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FIG. 15. Quotient of values derived for the bending mod- 
ulus using the model of [10] and that presented herein. 

Note that we have arrived at this expression by ignor- 
ing the transient component in the cxacit expression. This 
is valid only in the range where the polynomial (asymp- 
totic) solution dominates over the transients, a condition 
which is violated above p — 4, Xp/L — p / [p -\- 1) c:i 0.8. 

We expect then, for results obtained with the model 
in [10] systematically to underestimate the bending mod- 
ulus A by a factor ^1/2. In fact, wc can fit the published 
values of T(l) via nonlinear ^ fitting to arrive at a value 
for A. Unfortunately, the incorporation of the effects due 
to the left end of the polymer introduces a dependence on 
L, data for which are unpublished. Unable to elicit a re- 
sponse from the authors regarding the total lengths used 
in the experiment, we could only fit for both parameters. 
The analysis produced lower values than the model 
of [10]; however, is extremely wcakly-dcpendcnt on 
L and not a particularly advantageous method without 
this datum. In each case, the fit values of A were approx- 
imately twice that reported in [10]. 



VIII. CONCLUSIONS 

We have attempted to show that a systematic treat- 
ment of linearized clastohydrodynamics for filamentous 
biopolymers can be formulated with fruitful results. 
Specifically, we anticipate these methods should be use- 
ful in the design and analysis of dynamic experiments. 
Further, we have proposed a now technique (EHDIl) 
which exploits viscous hydrodynamics to extend the 
range of mechanical experiments of bending moduli to 
more flexible polymers. We expect this experiment to 
produce more accurate results when repeated with lower- 
amplitude driving, and thereby help determine conclu- 
sively the existence or nonexistence of scale-dependent or 
time-dependent clastic behaviors in biopolymers as well 
as the value of A in in vitro assays. It is our hope that the 
analysis associated with this experiment will also encour- 
age renewed interest in the problems of flagellar motion 
and slender-body hydrodynamics in general. 

Moreover, we have seen that attentiveness to equations 
of motion and boundary conditions for the elastica has 
measurable consequences, and that construction of the 
appropriate function space associated with these condi- 
tions leads to a pleasant union of mathematical and phys- 
ical consequences, quantifying our intuitions and relating 
transcendental equations to physical effects and experi- 
mental observables. We believe that the signiflcance of 
boundary conditions and the natural function space for 
the elastica has been overlooked in existing treatments of 
the dynamics and statistical mechanics relevant to exper- 
iments being performed and discussed by the community. 

Wc look forward to the extension of this analysis to ar- 
bitrary geometries, reflecting distortions beyond small or- 
der and hydrodynamics beyond the lowest-order approxi- 
mation of slender-body flow. Though the exact equation 
of motion is nonlinear, the Stokesian dynamic and the 
vanishing of successive derivatives remain, and wc ex- 
pect the mechanisms and effects which we have outlined 
here to persist. 

Natural extensions of this research include nonplanar 
geometries and the incorporation of twist. These would 
be complementary to recent work on the Hamiltonian dy- 
namics of twisted elastic rods [14, 25]. The dynamics of 
twist are especially intriguing in the light of recent work 
on twist-bend coupling [15] and the proposal that this 
coupling creates scale-dependent elasticity in actin [16]. 
Another intriguing experimental realization of nonplanar 
viscous clastohydrodynamics with twist is the supercoil- 
ing of fibers of mutants of the common bacterium Bacil- 
lus subtilis [24]. We are currently formulating the analysis 
appropriate to these promising applications, armed with 
the lessons learned in this investigation. 
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APPENDIX A: EIGENFUNCTIONS OF dt, 

If we equate functions related by the reflection a 
1 — a, there are 4(4 + l)/2 = 10 distinct eigenfunctions 
of 9^, determined by boundary conditions, for which this 
operator is sclf-adjoint. Each has an associated solvabil- 
ity condition for the eigenvalues k. We list the solvability 
conditions and the unnormalized eigenfunctions, indexed 
according to the conditions at the ends: 



(/) free : h^a: = Kxx = 
(c) clamped : h = = 
(h) hinged : h = h^x — 



{t) torqucd : 



hr, 







(Al) 

(A2) 
(A3) 
(A4) 



The general solution is Wk = ai cos ka + 02 sin ka + 
03 coshfca + 04 sinh/ca, where a G (0,1). The letters 
f,c,h,t denote the boimdary conditions at the left and 
right. Note that for special cases, the calculated Wk are 
merely a Fourier basis. 

• f - f : (A5) 
cos A; cosh A; = 1; 

Wfe = (sin ka + sinh ka) (sin k + sinh k) 
+ (cos ka + cosh ka) (cos k — cosh k) 
•c - c : (A6) 
cos fc cosh = 1; 

Wfe = (sin ka — sinh ka) (sin k + sinh k) 
+ (+ cos ka — cosh ka) (cos k — cosh k) 
•f - c : (A7) 
cosfccosh/c = —1; 

Wfe = (sin ka + sinh ka) (sin k — sinh k) 
+ (cos ka + cosh ka) (cos k + cosh k) 
•h - h : (A8) 
sin A: = 0; 
Wfe = sin ka 

•t - t : (A9) 
sin = 0; 
Wfe = cos ka 

•h - t : (AlO) 



cos fc = 0; 
Wk = sin ka 

•f-h: (All) 
tan A; = tanh A;; 

Wfe = (sin ka + sinh ka) (— cos A; + cosh A;) 
+ (cos ka + cosh ka) (sin A; — sinh A;) 
•f - t : (A12) 
tan A = — tanh A; 

Wfe = (sin ka + sinh ka) (sin A; + sinh k) 
+ (cos ka + cosh ka) (cos k — cosh A;) 
•c - t : (A13) 

tanh A = — tanh A; 

Wfe = (sin ka — sinh ka) (sin k -\- sinh k) 
+ (cos ka — cosh ka) (cos k — cosh A;) 
•c - h : (A14) 

tan A = tanh A; 

Wfe — (sin ka — sinh ka) (cos k + cosh A) 
-|- (— cos ka + cosh ka) (sin k + sinh A:) 



APPENDIX B: PROPERTIES OF THE Wk BASIS 

We define Wfe by the operator of which it is an eigen- 
function and by its boundary conditions: 



H = + boundary conditions. 



(Bl) 



Next we establish the self-adjointness of H, which we may 

express in Dirac notation as (Wfe|WWj) = (WWfe|Wj). 
This follows from the integration 

f'daWkdt (Wj) = Cdadt (Wfe) Wj (B2) 
Jo Jo 

+mdlw,\l-diWkdlw,\l 
+dlWkdiw,\l-dlWkWj\l 



J 

Jo 



dadt (Wfe) W 



3 ' 



□ 



where all the surface terms are seen to disappear for any 
choice of boundary conditions from appendix A. 

We can also establish the orthogonality of the members 
of a family since 



daWjWk 



' (^aWj^^ = k"" f dad^WjWk 
k Jq 



©1' 



' A* 
daWjWk . 



(B3) 



Thus /p daWjWk = for j ^ k. 

We justify considering only real, positive values of fc^ 
by observing 
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i'. 

Jo 



I daWkWl = J dad^WkWh 

+ [ dadlWkdlWl (B4) 
Jo 



so 



(B5) 



the quotient of two real, positive quantities, which is nec- 
essarily real and positive. 



APPENDIX C: EXACT SOLUTION TO EHD 
PROBLEM II 



the determinant of an 8 by 8 matrix to 0, an equation 
which can be written; 

Q = V{p,q) = 

— sin (gr) cos (p) cosh {q) cosh {p) 

— cos (g) sin (p) cosh {q) cosh (p) 
+ cos {q) cos {p) sinh [q) cosh {p) 
+ cos (g) cos (p) cosh (g) sinh (p) 

— sin (g) cosh (g) + cos (g) sinh (g) 

+ sin (p) cosh (p) — cos (p) sinh (p) (Dl) 

The ratio Xp/ {L — Xp) = q/p = p, set by the geometry of 
the experiment, fixes a diagonal line passing throngli the 
origin and intersecting the set of curves defined by this 
equation to define all the allowed k values, as illustrated 
in Fig. 13. 



The exact solution to EHDII can be written in some- 
what compact form at the cost of introducing new defi- 
nitions. Employing the parameter a = x/L, the rescaled 
length £ = L/£, and the constant zq = exp{—in/8) as in 
Eq. 5.11, we introduce 



and then write 

h{a, C) = hoo (a) + h\,{a, C) 



and 



(Cl) 



(C2) 



where /loo is the semiinfinite solution, and — > as 
£ — > 00. Explicitly, 

where 



(C5) 



APPENDIX D: SOLVABILITY CONDITION FOR 
DISCONTINUOUS EHDI 

Associated with the discontinuous EHDI problem is 
the set of eigenvalues of appropriate to the boundary 
and matching conditions. Like the statement sm{kL) = 
for determining the allowed Fourier k values for a doubly- 
hinged filament of length L, there is a solvability condi- 
tion for p = kxp , q = k{L — Xp) = kl, derived by setting 
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